Absolute emission rates of Spontaneous Parametric Down Conversion into single 

transverse Gaussian modes 
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We provide an estimate on the absolute values of the emission rate of photon pairs produced by 
spontaneous parametric down conversion in a bulk crystal when all interacting fields are in single 
transverse Gaussian modes. Both collinear and non-collinear configurations are covered, and we 
arrive at a fully analytical expression for the collinear case. Our results agree reasonably well with 
values found in typical experiments, which allows this model to be used for understanding the 
dependency on the relevant experimental parameters. 



INTRODUCTION 

For the last two decades, spontaneous parametric 
down-conversion (SPDC) has been the workhorse process 
for the generation of correlated photon pairs. These can 
easily be cast into maximally entangled states, which are 
useful e.g. for violating Bell inequalities [![, or for the 
investigation of other fundamental aspects of quantum 
mechanics. Progress in SPDC-based photon pair sources 
has allowed the emitted pairs from such sources to play 
a leading role in demonstration of quantum information 
techniques, and made its way to almost practical appli- 
cations like quantum key distribution. 

A number of studies have established the basic un- 
derstanding of SPDC @, 0, HI , based on energy conser- 
vation, and momentum conservation when participating 
light fields are treated as plane waves Q. 

Many of the more recent applications often necessi- 
tate significant manipulation and transport of the photon 
pairs; this can be achieved in a convenient way by guid- 
ing the light in single mode optical fibers. The basic idea 
of modeling SPDC in this regime is to map the optical 
modes propagating in the fibers into freely propagating 
modes of the electromagnetic field in the nonlinear con- 
version material, where they interact with a pump field. 
These freely propagating spatial modes can be described 
in good approximation by paraxial Gaussian beams, and 
any optimization strategy will involve some sort of mode 
matching of such interacting beams. 

Previous studies of SPDC light coupled into single 
mode fibers have focused in optimizing the coupling effi- 
ciency, defined as the ratio of photon pairs to single pho- 
tons that are observed because this is a quantity which 
can be measured easily in an experiment. This quantity 
is important for developing loophole-free tests of Bell's 
inequality [f| , heralded single photon sources 0, H, 0] , or 
simply sources of high pair brightness [Io| . 

So far, theoretical work in this area has focused mostly 
in such secondary parameters, and no closed expression 
for the absolute rate of photon pairs was available for typ- 



ical experimental configurations. This made it difficult to 
estimate whether a particular experimental source imple- 
mentation could be improved with respect to a particular 
figure-of- merit, be it total rate or spectral brightness. 

In this paper, we try to address this problem and de- 
rive an expression for the absolute rate of SPDC emis- 
sion from a bulk crystal into Gaussian modes. The work 
connects to earlier investigations of absolute SPDC rates 
with beams of finite diameter by Kleinman and Klyshko 
d, 0] . It was found there that the overall rate of pair pro- 
duction is independent of the spot-size of the pump beam 
[3] , and that the conversion efficiency of pump photons 
into correlated pairs integrated over all emission direc- 
tions [4] is in the order of 10~ 8 per mm for a typical 
non-linear material. The restriction on specific spatial 
modes defined by single mode optical fibers in the more 
recent applications, however, made it difficult to relate 
their results directly to experiments. Our description ap- 
plies both to Type I and II phase matching conditions, 
and covers collinear and non-collinear geometries impor- 
tant for the generation of polarization-entangled photon 



pairs 
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MODEL 

The basic process of SPDC can be understood as 
the spontaneous decay of a photon from a pump field 
into two daughter photons propagating in two — possibly 
different — target modes, with the process being medi- 
ated by a material with a nonlinear optical susceptibility. 
The physical implementation of SPDC utilizes the low- 
est order of the nonlinear susceptibility tensor m an 
appropriate material. While energy conservation allows 
the decay process to take place in many target modes, 
phase matching requirements need to be engineered to 
allow conversion to take place into any particular pair of 
directions. 

The physical model of the three interacting optical 
modes is depicted in figure [TJ We treat the pump beam 
and the target modes for the down-converted light as 
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FIG. 1: (Color online) Schematic of the down conversion 
model considered in this paper. Pump, signal and idler 
beams are treated as paraxial beams with a Gaussian trans- 
verse mode. The x-axis is perpendicular to the plane of 
the diagram. Coordinate systems of the signal and idler 



(a 



are tilted by angles 6 a ,9i with respect to the 



coordinate system of the pump. 

propagating paraxial beams with a Gaussian transverse 
profile. The beams overlap within a nonlinear optical 
crystal of finite length I, with surfaces normal to the 
propagation direction of the pump beam. Pump and 
target modes propagate in one plane, but need not to 
be parallel. We further assume that the three interacting 
modes overlap in a region without a significant variation 
of the transverse profile along their respective propaga- 
tion directions. This is a reasonable assumption for typ- 
ical Gaussian beam parameters and conversion crystal 
lengths used in experiments (Toj . 

We follow the tradition of referring to the target modes 
as signal (index s) and idler (index i) and choose coor- 
dinate systems where the 2 s ,i,p directions are parallel to 
the main propagation direction for each mode s,i,p (p 
refers to the pump mode) . The spatial mode function of 
the electrical field for each of the modes can be written 
as 



U(x,y) = e 



ikz 



(1) 



where k denotes the z-component of the corresponding 
wave vector, W the Gaussian beam waist parameter, 
and x,y 7 z refer to the corresponding coordinate system 
for each mode. To simplify the overlap calculations, we 
use normalization constants a for the envelope functions 
U(x, y) such that 



a 2 / dxdy \U(x, y)\ 2 = 1 



(2) 



in their corresponding coordinate systems, which implies 

/ 9 

(3) 



irW, 



p.s.i 



We note that the spatial mode function g(r) fulfills 
Maxwell's equations only approximately. For the calcu- 



lations presented below, however, this poses no problem. 
Furthermore, the dispersion relation connected with this 
mode function has the confinement correction: 
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(4) 



Again, for practical beam diameters W of about 100 
wavelengths considered in this paper, this correction term 
is small enough to safely neglect it. 



Pump mode 

The pump mode is aligned with the main coordinate 
system x, y, z, and treated as a classical monochromatic 
field of amplitude E p . We further assume (as is custom- 
ary) that we have no significant depletion of the pump 
in the downconversion crystal. The electrical field of the 
pump can thus be written as 



E p (r,t) = MEW(r,t) + EH(r,t) 



(5) 



[E°e p g p (r)e- 



+ c.c. 



with a polarization vector e p , and a corresponding an- 
gular frequency lo p . Using the normalization expression 
@ , we can connect the electrical field amplitude E p with 
the optical power P in the pump beam: 



\EP\ % =c?-™- 
1 Pl ^eon p c 



(6) 



with the refractive index n p for the pump field, the elec- 
trical field constant eo and the speed of light c in vacuum. 



Target modes 

First, we need to take care of the possibly non-collinear 
propagation of the target modes with respect to the 
pump. By introducing target mode angles # Sj i, and us- 
ing an orientation as indicated in figure [TJ we express the 
spatial coordinates of the target modes in terms of the 
main coordinates x, y, z: 









cos6L 




(7) 



To arrive at a rate of photon pairs generated via SPDC, 
we use the quantized field operators. We do that by in- 
troducing a quantization length L in propagation direc- 
tion for clarity in the counting of modes, and postulate 
periodic boundary conditions; later we will drop this re- 
quirement. Following the notation in equation ([5]), the 
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electrical field operators take the form 
E Sti = i[feW(r,t)+fey(r,t)] 



kTl V n ^ e ° 

+ /l.C. 



(8) 



Here, e s .j indicate the polarization vectors, and n Sj i and 
w Sj j the corresponding refractive indices and angular fre- 
quencies of the target modes. The modes are indexed by 
the scalar moduli fc Sj i, and the corresponding full wave 
vectors in the pump coordinates are given by 



fc S) i(=Fsin0 S) iej, + cos9 Sti e z ) 



(9) 



The longitudinal wave vector components serve as a 
complete discrete mode index k St i — 27tm s .i/L with in- 
teger numbers vn s ^. The coefficients before the raising 
and lowering operators in cq. (jHJ) are chosen such that 
the free field Hamiltonian Hq for the target modes takes 
the usual form 



"k s . 



■ a k s 



(10) 



Interaction Hamiltonian 

The SPDC process is enabled by a nonlinear optical 
material whose presence is described by the Hamiltonian 
Hi, written in the interaction picture with time depen- 
dence of the raising and lowering operators [l3[ : 



(2) 
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d J dxdy J dz 

-oo -1/2 fes ' fci 



x e- iA »\(r)g:(r)g*(r)rtSt)al(t) + h.c.(ll) 

We have assumed a crystal of infinite transverse (x, y) 
extent, which is justified when the beam diameters are 
much smaller than the crystal dimensions. We introduce 
a frequency mismatch Auj = lu p — u> s — uji . The effective 
non-linearity d captures the contraction of the nonlinear 
susceptibility tensor with the corresponding polarization 
vectors {2d = e p x^ : e s ej) (la ]. With this notation, 
the type of phase matching condition (type I or II) is 
reflected in an appropriate effective non-linearity d. 

Most of the scaling aspects of the parametric down 
conversion process connected with the geometry of the 
interaction are determined by the overlap integral $(Afc) 



of the three mode functions g p . s ^(r) in the crystal: 
$(Ak) = Jdzjdydxg p (r)g*(r)g*(r) 

fdz fdydxe^Upi^Usirp^r). (12) 



In this expression, Ak = k p — k s — describes the wave 
vector mismatch. Since pump and target modes are de- 
fined in the y-z plane, there are no wave vector compo- 
nents in the ^-direction and hence Ak x = 0. Carrying 
out the integration in the transverse directions (x, y) we 
arrive at 
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(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 



The exponential term before the residual integral in 
cq. (|13p represents the approximate transverse wave vec- 
tor mismatch. This term can be ignored only if one of 
the beams is infinitely large {W p , s ,i oo), or if there is 
perfect transverse phase matching. 

The residual integral along z in (|13p can be re-written 
in a form that allows for a physical interpretation. Wc 
introduce where 



1/2 

& z := J dze 

-1/2 



Hz 2 +izK 



= I- jdue au cos(Atpu). 
" 



(20) 



(21) 



The phase mismatch is now defined as A(p := Kl/2. 
The argument 5 := \/~Hl/2 in the exponential can be 
viewed as a "walk-off" parameter due to non-collinear 
mode propagation. This parameter is useful for identi- 
fying a thin crystal and a thick crystal regime [6(. In 
our model, these regimes refer to the physical boundary 
conditions imposed on the interaction volume by the ge- 
ometry of the pump and target modes. 
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FIG. 2: (Color online) Longitudinal overlap function $ Z /Z as 
a function of the total phase mismatch Atp = Kl/2 over the 
crystal for various walk-off parameters H. For H — 0, the 
typical sine shaped spectral distribution is revealed, whereas 
for large walk-off parameters 3 > 1 the phase matching con- 
dition is determined by the overlap region formed by pump- 
and target modes, and develops into a Gaussian distribution. 



In the thick crystal regime with a large walk-off param- 
eter (3 > 1), the overlap integral & z depends mostly on 
the characteristic beam parameters W PjS> i and not much 
on the physical boundaries of the non-linear material. 
For S — > oo the length of the crystal ceases to play a role 
altogether: 
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l h Erf(s) 

^E.riE) 



(22) 



The thin crystal regime refers to a small walk-off pa- 
rameter, 2 C 1, so that the characteristic beam pa- 
rameters have almost no influence on <I> Z . In particu- 
lar, this applies for collinear arrangement of all modes 
(9i = S = 0), where 5 = 0. In this case, K = Ak z , and 



$ z = I sinc(A<£>) 



(23) 



Spectral Emission Rate 

To obtain absolute emission rates, we make use of 
Fermi's Golden Rule as an expression for the transi- 
tion rate R(k s ) between the initial vacuum state \i) = 
|0fc s ,0fc ; ), and a final state |/) = a\ 10^,0^) with the 
mode pair k s , ki populated with one photon each. Fermi's 
rule applies for asymptotic scattering rates, so the rela- 
tion between k s and ki is fixed by energy conservation: 



Ato 



k a — kj — — 



(24) 



We first evaluate the transition rate R(k s ) to a fixed 
signal target mode k s . The density of states p per unit of 
energy AE — TiAuj is extracted out of a quasi-continuum 
of states for the mode ki\ 



p{AE) = 



A? 



Ah d(hAw) 



L rij 
2ir he ' 



(25) 



where Am/Aki — L/2tt denotes the number of modes 
per unit of wave vector component ki . 

With the transition matrix element expressed in terms 
of the overlap integral $(Ak), 



(f\H I \i) = d^^^E^(A^ : 



n s m L 
the transition rate is then given by 
2tt 

\d a s aiE°${Ak) 



R(k s 



(f\Hi\i) P (AE) 



U! s Uli 

nlnicL 



(26) 

(27) 
(28) 



The spectral emission rate per unit of angular fre- 
quency Lo s is obtained by multiplying R(k s ) with the 
number of modes k s in a unit interval of ui s , which is 
Ln s /2irc. We finally arrive at 



dR(ui s ) 
duo* 



d a sai £°$(Ak) 



lixn^rii 



(29) 



At this point, the earlier introduced quantization 
length L has vanished as expected. 



This reveals the well-known influence of the longitu- 
dinal phase mismatch on the downconversion spectral 
properties 0] . Figure [5] shows the overlap contribution 
<& z /l as a function of the phase mismatch A<^ for vari- 
ous walk-off parameters S. We note that as S becomes 
large, the spectrum becomes Gaussian-like. If 3 is iden- 
tified as the degree of (non)collinearity, it suggests that 
the SPDC spectral profile far from collinear emission can 
be very broad. In other words, for the same beam pa- 
rameters W S) i ]P , the bandwidth of collinear emission will 
be narrower than for non-collinear emission. 



Total Emission Rate 

We now can determine the total pair generation rate 
by integrating the spectral rate density over all frequen- 
cies ui s . Assuming that the overlap $(Ak) is only non- 
vanishing over a small range of frequencies ui s , the total 
pair generation rate can be written as 



Rt = 



dctMiE® 



2irn„rii 



cM$(Ak)| 2 



(30) 
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The dependency of <I>(Ak) on uj s can be quite involved, 
as in the non-collinear case (9i tS ^ 0) both Ak y (L0 s ) and 
Ak z ((jj s ) must be considered. However, the alignment 
criteria for most experimental setups assume perfect lon- 
gitudinal phase matching to arrive at the collection an- 
gles for degenerate downconversion 2, 11 1 . In these col- 
lection directions, the target mode angles 6 Sj i are equal. 
Furthermore, the typical experiments use identical col- 
lection mode diameters for signal and idler, W s = Wi 
[icj |. Under these two conditions, the phase mismatch 
Aip depends only on Ak z . 

We now consider the exponential term for the over- 



lap $ in eq. (|T3| that contains Ak y . For experiments 
where light centered on the degenerate wavelengths with 
a small bandwidth is collected (w 2 nm on either side 
of the center [l(|), we will assume perfect transverse 
phase matching. A complete treatment with non-zero 
transverse phase mismatch requires a numerical proce- 
dure 



1|. 



With perfect transverse phase matching, we can carry 
out the integration in eq. (|30|) by re-parameterizing the 
frequencies of the signal and idler about the degenerate 
SPDC frequency: uj s — ^ — 5^ and = ^j- + S^. We 

approximate uj s uji w -f- by ignoring terms 0(6^). From 
energy conservation eq. (|24|) and the phase matching con- 
dition 
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FIG. 3: Variation of the spectral integral S with walk-off 
parameter S. Since the absolute emission rate is proportional 
to S, the largest absolute rate is obtained in the thin-crystal 
limit with H = 0. 



limit, closed form expressions for the spectral and total 



rates are 

dR(u s ) 



du> s 



!?-/ 



2d 2 w 2 PPsinc 2 (A/c z Z/2) 
97rn p n s nie c 3 (ttW 2 ) 
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and 



(35) 
(36) 



Ak z = n s oj s cos 9 S + m&i cos 9i 



(31) 



we obtain a dispersion relation between du s and d{Ak z ): 



d(Ak z ) 



(n-i cos 9i — n s cos S ) 



dLUe 



(32) 



The emission rate can now be integrated over the lon- 
gitudinal wave vector mismatch Ak z , 



Rt — 



d 2 {a l a s E° p ) 2 u J 2 p 



Acn s rii (2-7r) (n^ cos 9i — n s cos 9 S ) 



\t>(Ak)\ 2 d(Ak z ). 

(33) 

Effectively, this is the pair emission rate for all allowed 
wavelengths in the direction defined by our paraxial 
beams. If we recall that the pump has a Gaussian en- 
velope, and choose all beam characteristics to be equal 
(W p = W s = Wi), the rate Rt finally can be written as 



Rt 



Ad 2 Plujl 



37m p n s n i eoc 2 (7rW / ?)(l + cos# 4 2 + cos6> 2 ) 



S. 



(m cos 9i — n s cos 9 S ) 



(34) 



with the spectral integral S := / — ^ — — d(Ak z l/2). 
The absolute emission rate is proportional to S, which 
has a dependence on the value of the walk-off parameter 
S as shown in figure [3] The spectral integral assumes 
its largest value S — tt in the thin crystal limit. In this 



Dependence of Emission Rate on Beam Waists 

Although it is convenient to set all beam waists to be 
equal, this is not necessary. In fact, it can be shown that 
this choice does not maximize the total emission rate 
for a given optical pump power. Carrying out the more 
general derivation to arrive at an expression similar to 
(|5rJ|) . the dependency on the various beam waists W p , W s 
and Wi (again in the thin crystal limit) can be written 
as 



Rt oc 



W 2 W 2 W 2 { 



1 

H'.T 



1 

w 2 



i 



(37) 



To develop an alignment strategy, we may assume that 
the collection modes are identical (W s — Wi — W), but 
we re-express the pump waist as W p = jW, so we obtain 



Rt oc Rt oc 



1 



IF 2 (i + 2 7 ) 2 



(38) 



This relationship is illustrated in figure IU and exhibits 
a maximum of Rt for 7 = -4=. For 7 = 1, the emission 
rate is about 12% lower than the maximum value. This 
suggests that experimental setups that are designed with 
equal beam waists for pump and collection modes may 
be further optimized, and the simple argument of max- 
imizing a mode overlap (loj with matching beam waists 
does not hold. 
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FIG. 4: Dependence of the total pair rate Rt on the ratio 7 
between target and pump waist. The maximum emission rate 
can be expected at 7 = 1/ v2- 



PHYSICAL INTERPRETATION AND 
COMPARISON TO EXPERIMENTS 

While Rt is proportional to crystal length I, the spec- 
tral rate dR{uj s )/duj s is proportional to the square of a 
sine function. This is in agreement with results from 
previous work [3J. However, our expression reveals de- 
pendencies on other factors, namely emission geometry 
and pump spot-size. The expression for Rt reveals that 
the emission rate is higher in a collinear geometry com- 
pared to a non-collinear case. This can be intuitively 
understood because the non-collinear case has a smaller 
interaction volume. 

We find that the absolute rate Rt is proportional to 
the square of pump frequency since we re-parameterize 
the signal and idler about the degenerate frequency, so 
downconversion efficiency can be improved with shorter 
wavelength pumps. Both spectral and total emission 
rates are inversely proportional to the mode area of the 
beams, in contrast to previous studies which showed that 
the total SPDC cannot be enhanced by focusing [3L flj|. 
There, however, SPDC emission was not considered for 
specific transverse modes. The dependence of emission 
rates on the mode area have been reported in a previous 
analysis of SPDC in waveguide structures This is 

not surprising because the emission into paraxial beams 
is essentially the same problem as SPDC in waveguides, 
where the target modes are quantized in one dimension 
only. We note that our equation (|35p is similar to the 
equation obtained in reference [l6j |. 

We should not draw the conclusion, however, that 
SPDC emission into single transverse modes can be ar- 
bitrarily enhanced by tight focusing. Our model is only 
valid in cases where the transverse profile of the beams 
does not vary significantly over the crystal length. For 
an optimization study of focus size on SPDC emission we 
refer the reader to [9j . 

For explicit comparison of eq. (|34[) with experimental 



values, we consider our experimental setup (similar to 
that used in reference [13] )• I n this experiment a pump 
beam (beam waist, W p = 82/im) at a wavelength of 351.1 
nm is incident on a 2 mm thick /3-Barium Borate (BBO) 
crystal. Two single mode fibers are used to collect de- 
generate downconverted photons, which is estimated to 
have an external emission angle of 3.1°. The collection 
modes also have beam waists of W s ,i — 82/xm. 

Specifically for BBO the effective non-linearity is given 
by d = c^22 cos 2 9 C cos 30 c . The angle between pump wave 
vector and crystal optical axis is 9 C — 49.7°, while the az- 
imuthal angle is <fi c — 60°, resulting in an effective non- 
linearity of 9 x 1CT 13 m/V (d 22 = 2.11 x 1CT 12 m" 1 V" 1 
according to [HI). The observed pair rate is approxi- 
mately 800 pairs mW^V 1 with a pair-to-singles ratio 
of 0.23. 

The walk-off parameter for this setup is S = 0.933, 
indicating that the overlap integral is intermediate be- 
tween the thin and thick crystal limits. The largest ob- 
servable rate according to our model is 2(0.23 x Rt) = 
1100 mW^V 1 . The additional factor of 2 is used be- 
cause in experiments, the geometry is used to collect 
downconversion emission in two decay paths. 

The source of the discrepancy between experiment and 
our model is hard to identify. The assumptions used in 
the model make it an overestimate, primarily in the re- 
parameterizing of signal and idler frequencies about the 
degenerate wavelength. Experimentally, there are several 
sources of uncertainty, the main one being the difficulty 
in establishing pump power very accurately. For exam- 
ple, the average observed value was arrived by measuring 
the power using two different power-meters (a Newport 
Model 818-UV reported 11.7 mW while a Coherent Field- 
master reported 9 mW) . The uncertainty in pump power 
estimation, however, is not sufficient to make the ob- 
served result compatible with the calculated value. 

According to the model, the conversion efficiency into 
Gaussian transverse modes for our experimental setting 
will be 3 x 10~ 12 mm -1 of crystal length. Other exper- 
imentally reported rates in the literature reveal similar 
downconversion efficiencies [H, [13, EH ■ The total conver- 
sion efficiency of SPDC was found to be on the order of 
3xl0~ 8 mm -1 sr _1 by Klyshko [1] (for degenerate SPDC 
with a 500 nm pump wavelength). Experimentally, our 
collection angle is 3.3 x 10~ 5 sr. If we convert our units 
to be comparable with Klyshko's result, we obtain an 
efficiency of 7 x 10~ 8 mm -1 sr -1 . 

CONCLUSION 

In conclusion, we have presented expressions that pro- 
vide absolute values for the rate and bandwidth of cor- 
related pairs emitted in bulk crystal SPDC which are in 
a single Gaussian mode. These modes may be defined 
by the collection profile of single mode fibers, selecting 
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emission in a specific pair of directions. The single mode 
treatment reduces the complexity in the final expression 
for the rate equations. 

We find that the expression for absolute rates given by 
the model are slightly larger than experimental observa- 
tions. The model may thus be treated as an idealized 
case for the total pair emission rate. The small differ- 
ence between experimentally observed rates and predic- 
tions according to the closed expression for Rt, however, 
suggests that experimental setups using single mode col- 
lection fibers (e.g. HEH) operate close to the optimal 
limit. 

Substantial increase of the emission rates are to be ex- 
pected from larger non-linearities, since emission rates 
are proportional to d 2 . Small mode diameters are also 
expected to enhance emission rates, as has been con- 
vincingly reported for SPDC experiments using waveg- 
uide structures I 
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|22j, and a similar theoretical 
analysis [16| for those cases. Overall spectral brightness 
will be improved by combining larger non-linearities with 
collinear mode confinement in longer structures. Even 
then, however, the spectral width is still ultimately deter- 
mined by the longitudinal wave vector mismatch. This 
indicates that very dramatic improvements (by several 
orders of magnitude) to the generated pair rate in a nar- 
row bandwidth necessary for addressing atomic systems 
is not very likely to be expected from bulk crystal emis- 
sion. 
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